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[57] ABSTRACT 

Doppler in a communication system operating with a 
multiple differential phase-shift-keyed format 
(“MDPSK”) creates an adverse phase shift in an incom- 
ing signal. An open loop frequency estimation is de- 
rived from a Doppler-contaminated incoming signal. 
Based upon the recognition that, whereas the change in 
phase of the received signal over a full symbol contains 
both the differentially encoded data and the Doppler 
induced phase shift, the same change in phase over half 
a symbol (within a given symbol interval) contains only 
the doppler induced phase shift, and the Doppler effect 
can be estimated and removed from the incoming sig- 
nal. Doppler correction occurs prior to the receiver’s 
Final output of decoded data. A multiphase system can 
operate with two samplings per symbol interval at no 
penalty in signal-to-noise ratio provided that an ideal 
low pass pre-detection filter is employed, and two sam- 
ples, at J and $ of the symbol interval T s , are taken and 
summed together prior to incoming signal data detec- 
tion. 

31 Claims, 13 Drawing Sheets 
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DOPPLER-CORRECTED DIFFERENTIAL 

DETECTION SYSTEM 

ORIGIN OF THE INVENTION 5 

The invention described herein was made in the per- 
formance of work under a NASA contract, and is sub- 
ject to the provisions of Public Law 96-517 (35 USC 
202) in which the Contractor has elected not to retain 
title. 10 

TECHNICAL FIELD 

This invention relates to an improved communication 
system which is provided with a Doppler-corrected 
differential detector located in the communication sys- 15 
tern’s receiver. 

BACKGROUND ART 

Digital data communication systems employing a 
carrier that is modulated by data at a transmitter and 20 
demodulated at a receiver to decode the transmitted 
data are well known. Several different types of modula- 
tion and demodulation are known and used. This inven- 
tion pertains to a type of modulation-demodulation 
which employs data transmission through multiple dif- 25 
ferential phase-shift-keying (“MDPSK”) at the modula- 
tor and a comparable differential detection at the sys- 
tem’s receiver or demodulator. 

Modems, a contraction for modulators and demodu- 
lators, are commonly used at both ends of a data com- 30 
munication transmission system. When one modem 
physically moves relative to another modem in a trans- 
mission link, the transmitted signals experience a Dop- 
pler effect. Typically Doppler effects are experienced in 
an earth to space satellite link and cause an effective 35 
change of frequency of a received signal due to the 
relative velocity of the link’s transmitter with respect to 
the link’s receiver. The Doppler effect must be rectified 
in order to achieve efficient data transmission and re- 
ception. 40 

A prior art search related to the invention discloses 
some typical approaches to Doppler correction. Not 
one of the references located in the search, however, 
discloses or suggests an open-loop Doppler frequency 
estimation derived from an incoming signal. Moreover, 45 
open-loop correction for a Doppler effect in a MDPSK 
system having an improved filtering and differential 
detection configuration is not disclosed in the prior art. 

The following search-located U.S. Pat. Nos. include: 
Heard et al 4,520,493, Poklemba 4,419,759, Poklemba et 50 
al 4,472,817 and Mehrgardt et al 4,663,595. 

Heard ’493 discloses a complex closed-loop Doppler 
correction scheme that decodes and then attempts to 
correct for the next expected Doppler effect. Such 
closed loop tracking is not feasible in certain Doppler- 55 
related-applications e.g., a multipath fading environ- 
ment. Poklemba et al ’817 discloses open loops but does 
not disclose Doppler correction in a differential detec- 
tion system. Poklemba ’759 is a typical phase-locked 
loop system. Mehrgardt ’595 is of marginal interest only 60 
in that it discloses delays and summers in a form typi- 
cally known in the modulator/demodulator art. 

A multipath fading . environment creates uncertain- 
ties in the carrier phase and/or frequency. These uncer- 
tainties essentially eliminate the known techniques 65 
which use a phase and/or frequency tracking synchro- 
nization loop. Moreover, the above-noted uncertainties 
may be further masked by phase jitter on the carrier 


2 

reference. This additional transmission impairment adds 
further complexity to the system operation. Differential 
detection employs, as a demodulation reference, the 
phase and frequency of the carrier corresponding to the 
previous transmitted data symbol. The deleterious ef- 
fects in a multipath fading environment are largely 
avoided by that differential system approach because 
there is no requirement that the carrier must be syn- 
chronized during each present symbol interval. 

Another known transmission impairment is intersym- 
bol interference (ISI). It is known to use pulse shaping 
so that the ISI at a receiver is at a minimum at the sam- 
pling point for data detection by the receiver Splitting 
signal filtering between the transmitter and receiver 
locations is a known way of achieving pulse shaping 
and minimum ISI. 

For MDPSK systems it is common to strive for what 
is known as a 100% raised-cosine spectrum at the de- 
coder’s input. In an MDPSK system such shaping will 
achieve a satisfactory signal to noise ratio for center 
sampling at the decoder or demodulator. Various types 
of multilevel coding and decoding schemes have been 
introduced in efforts to achieve an optimum signal to 
noise ratio for such center sampling. Included within 
these multilevel schemes is the so-called duobinary 
signaling, a partial response signaling scheme as taught 
in A. Lender, “The Duobinary Technique for High 
Speed Data Transmission,” IEEE Transactions on 
Communication Electronics, Vol. 82, May 1963, pp. 
214-218. 

STATEMENT OF THE INVENTION 

The invention features an open loop frequency esti- 
mation derived from a Doppler-contaminated data- 
modulated input signal. A Doppler estimate is obtained 
in a doppler-estimating circuit and that estimate is fol- 
lowed by Doppler correction prior to the receiver’s 
decoding of the modulated data. An overall Doppler- 
corrected differential detection system of simple and 
efficient design is thus provided. 

One of the key features of the invention is based upon 
the inventor’s recognition that whereas the change in 
phase of the received signal over a full symbol contains 
the sum of the differentially encoded data and the Dop- 
pler induced phase shift, the same change in phase over 
half a symbol (within a given symbol interval) contains 
only the Doppler induced phase shift. Thus, by proper 
open-loop processing, including post-detection integra- 
tion (whose duration is related to the rate of change of 
the Doppler), the Doppler effect can be estimated and 
removed from the incoming signal. 

It has been further discovered that a multiphase sys- 
tem can operate with two samplings per symbol interval 
at no penalty in signal-to-noise ratio. This discovery is 
possible provided that an ideal low pass pre-detection 
filter is employed, and two samples, at \ and | of the 
symbol interval T* are taken and summed together 
prior to decoding. 

BRIEF DESCRIPTION OF THE DRAWING 

FIG. 1 depicts a prior art differentially coherent re- 
ceiver employing a full symbol delay; 

FIG. 2 is a basic block diagram of an open-loop Dop- 
pler-corrected differentially coherent receiver system in 
accordance with the invention; 

FIG. 3 depicts an I-F implementation of the invention 
of FIG. 2; 
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FIG. 4 depicts a baseband implementation of the 
invention of FIG. 2; 

FIG. 5 is a set of curves useful in demonstrating the 
performance of the invention; 

FIG. 6 is another set of curves also useful for demon- 5 
strating the invention’s performance; 

FIG. 7 is a simple building-block diagram depicting 
the Doppler-corrected receiver in complex form, in a 
tracking mode; 

FIG. 7A supplies a more detailed form of the Dop- 10 
pier estimator of FIG. 7; 

FIG. 7B depicts a fast acquisition circuit for Doppler 
correction; 

FIG. 7C depicts a fast acquisition circuit which is 
equivalent to that depicted in FIG. 7B; 15 

FIG. 8 is a block diagram, employing the building 
blocks of FIG. 7, 7A and 7C, to form a dual mode 
Doppler-corrected differential detector in complex 
form; 

FIG. 9 is the dual mode Doppler-corrected MDPSK 20 
receiver of the invention using quadrature components 
of a received baseband signal; 

FIG. 10 is an eye diagram that shows sample points 
for zero intersymbol interference for the system of FIG. 

9; 25 

FIG. 11 is a mathematical model which is useful in 
describing a more complex feature of the invention; 

FIG. 12 is a mathematical equivalent of FIG. 11; 

FIG. 13 is a graph depicting the improvement in 
results achieved by the invention in the absence of time 30 
jitter; and 

FIG. 14 is a graph similar to that of FIG. 13 depicting 
results in the presence of ten and twenty percent time 
jitter. 

DETAILED DESCRIPTION OF THE 35 

INVENTION 

FIG. 1 is a typical block diagram of a known differen- 
tial detector in a MDPSK receiver. At input terminal 5 
an input signal s(t) which may include modulated data 40 
and Gaussian noise is split into an in-phase (I) and a 
quadrature phase (Q) channel 10 and 20 respectively. 
The quadrature phase channel 20 includes a 90 degree 
phase shifter 15. 

Each channel includes a full symbol interval delay 45 
provided by two half-symbol delay networks 31, 32 and 
41, 42. Each network thus provides a one-half symbol 
delay time. Multiplying junctions 35, 45 multiply the 
undelayed and delayed signals together to form a basis 
for comparing two different time-related phase angles 50 
of the input signals applied to those multiplying junc- 
tions 35, 45. The output of each multiplying junction is 
applied to a detector 50, which detector determines the 
angular difference between s(t) and s(t) delayed by one 
symbol and applies a decoded output signal an output 55 
lead 6. Both in-phase (I) and quadrature phase (Q) chan- 
nels 10, 20 are required to represent the well-known 
phasor diagrams of the prior art. 

FIG. 2 includes the broad concept of the Doppler- 
corrected differential phase detection of the invention. 60 
A key feature of the invention is provided by a Doppler 
estimator circuit 100 which requires input signals from 
both the in-phase channel 10 and the quadrature chan- 
nel 20. The in-phase channel 10 supplies a pair of input 
signals to circuit 100 in the form of an incoming signal 65 
s(t) and a delayed version thereof. Doppler estimator 
circuit 100 receives the incoming signal s(t) over lead 90 
and also receives a one-half (T/2) delayed version of 
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that incoming signal over lead 101, which lead connects 
the output of the T/2 delay circuit 31 to Doppler esti- 
mator 100. Simultaneously the quadrature channel 20 
applies an input signal to circuit 100 in the form of a 
phase shifted s(t) and a delayed version thereof. Being 
delayed by less than a full symbol interval, the input 
signals from channels 10 and 20, into Doppler estimator 
100 have the same data-representing phase; but, with 
proper processing, the Doppler-induced change in 
phase in the incoming signal s(t) can be estimated in 
circuit 100. 

Multipliers 35 and 45, operating as before, are con- 
nected as in FIG. 1. Doppler removal is done in both 
the channels and it will now be described. 

Between the multipliers 35, 45 and decoder 50 are 
located Doppler removal circuits 125 and 150. Each 
Doppler removal circuit receives the estimated Dop- 
pler signal from estimator 100 and that estimated signal, 
prior to signal decoding in detector 50, is removed by 
Doppler removal circuits 125, 150 in the (I) and (Q) 
channels 10, 20, respectively. With the Doppler effect 
nullified, detector 50, in known fashion, delivers a data- 
representing output signal on lead 6 in accordance with 
standard MDPSK decoding techniques. 

The detector 50 may be of any suitable and known 
type. Such detectors include the type which employ a 
maximum-likelihood decoding algorithm. Further pro- 
cessing of the output signal on lead 6 depends upon 
whether the transmitted data is coded or uncoded. In 
the case of uncoded data, a decision on the output is 
conventional as known for FIG. 1. For coded transmis- 
sion, the output is used directly (in practice, it is Q-BIT 
quantized) as a soft decoding metric. 

A more detailed and rigorous description of the in- 
vention is now presented with reference to FIG. 3. In 
FIG. 3, the Doppler-contaminated input signal s(t) is 
shown at input terminal 5. That signal is defined as: 

VI (0 

s(0 = N 2P sin [(tuo + Ao>)r + 0(0] 

where &>0 is the carrier, and Aco the Doppler shift and 
0(t) is the data modulated phase information in the form 
of a predetermined angle shift and some arbitrary phase 
shift that may be neglected. The input signal will also 
contain Gaussian noise n(t) at lead 7, which noise is 
defined as: 

\ I — (2) 

n{t) = N 2 {«c(0 cos (too + Aa >)t -j- n^r) sin (a>o + A w)/} 

where P denotes the average signal power in watts, 
Aa) = 27r0f is the radian Doppler shift, 8(t) is MPSK 
modulation with symbol rate 1/T S , and nc(t), ns(t) are 
low-pass Gaussian noise processes with one-sided 
power spectral density No w/Hz. 

Input leads 5 and 7 are connected to summing junc- 
tion 14, which junction, in turn, is connected to a band 
pass IF filter 33 which functions to limit the noise in the 
summed signal. After passing the sum of s(t) and n(t) 
through the IF filter 33 (which for simplicity of analysis 
is assumed not to distort s(t», the resulting signal, y(t), 
is then processed in parallel by the in-phase (I) and 
quadrature phase (Q) channels 10 and 20. These two 
channels produce signals suitable for data detection and 
Doppler estimation in accordance with the present in- 
vention. 
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In the Doppler estimation circuit 100 of the receiver, 
the in-phase and quadrature input components (delayed 
Ts/2 by delay circuits 31, 41) are multiplied at multipli- 
ers 21, 22 with the undelayed input itself. After zonal 
filtering by IF filters 34, 36 the second harmonic of the 
carrier has been removed, leaving the following signal: 


xdO = yiW - Ts/2) 

= P cos [(coo + Aco)7V2 + 0(f) - 6(t - Ts/2)} + N^t) 

**<0 = - Ts/2) 

= Psm [(«o + Aui)Ts/2 + 0(0 - 0(r - Ts/2)] + N&) 


NA 0 - nMt - Ts/2) 4- n{t - 7^2)*/) 4- n{t)n(t - Ts/2) (4) 

Ns(t) = «(0^90(^ - Ts/2 ) 4- *9oO ~ V2W0 + «(0«9o(' - 7s/2) 

In the above equations (3) and (4) N c (t) and N^t) are 
noise terms, and the subscript 90 refers to a 90° phase 
shift in the signal having that subscript. Making the 
assumption, as is typical in differential detection analy- 
ses, that 0) o T s /2 is an integer multiple of 2 tt, then Equa- 
tion (3) may be simplified to 


xdt) = P cos [Aa>Ts/2 + 0(0 - 8(t - Ts/2)] + N^t) (5) 

x^0 = P sin [Aco7V2 4- 0(0 - 0(t - Ts/2)) + N£t) 

Since the MPSK modulation 0(t) is constant over a 
symbol time, then 0(t)— 0(t— Ts/2)=0 in the latter half 
of every symbol interval. 

In FIG. 3, switches 81, 82 denote time sampling by 
any conventionally well-known sampling circuits. If we 
sample x c (t) and x^t) at U=[(4k+3)/4]T 5 , when k is 
equal to 1,2, 3 etc., then these sample values denoted by 
x c k and Xjjfc, respectively, become 


-continued 

y s — sin Ao)T s 

gives a pair of quadrature signals that can be further 
processed and used to remove the Doppler content 
from the in-phase and quadrature data detection chan- 
nels, 10, 20 respectively. Before describing that re- 
moval, however, the data detection portion of the sys- 
tem will be described in further detail. 

In the data detection portion of the receiver, the 
in-phase and quadrature input components (delayed by 
T 5 via circuits 31, 32 and 41, 42) are correlated with the 
input itself to produce (after zonal filtering in circuits 
51, 52) 

*c'(0 = y(w - T s ) (9) 

= P cos [Aw r, + 6 ( 1 ) - 6 (t - T s )] + N c '(t ) 

*sV) = - t s ) 

= P sin [AojTs + 0(0 - 0(r - T s )] + N s '(t) 


N C V) = nMt - T s ) 4- n(t - T s )s(t) + n(t)n(t - T s ) (10) 

Ns( 0 = n(t)s 9 o(t - T s ) 4- ~ T s )s( t) + n(t)ngo(t - T s ) 

and, in accordance with out previous assumption, we 
have made use of the fact that Ao))T 5 is an integer multi- 
ple of 277. Switches 53, 54 are used to sample Xc(t) and 
x*(t) at U=[(2k— 1)/2]T S (i.e., the midpoints of the sym- 
bol intervals). See FIG. 10 for the location of midpoint 
sampling. 

We obtain corresponding noise free sample values of 


x c k = P cos Ao^Ts/2 + N c k (6) 

x 5 k = P sin Aa>7V2 4- N s k 

45 

The signals x c /tand x 5 £are sampled by switches 81, 82 
respectively and these samples are amplified by gain 
circuits 91, 92 each of which provide a gain term of 
1 — a. Equation 6 shows that the data term has disap- 
peared and we are left with the Doppler estimation 50 
signal. Finally, filtering each sample sequence to reduce 
undesired phase noise, an estimate Ao)T 5 of the normal- 
ized Doppler AofTjCan be obtained. The required filter- 
ing is performed using digital filters 93, 94 with gain 
coefficient a (0<a<l). In such a case, the Doppler 55 
estimate, Ao>Ts, is given by 


Xck' — P COS (Aoj7V 4- ABk) 4- N c k' 


x s k‘ — P sin (Ao )Ts 4- A 0k) + N s k' 


Aa}T s — 2 tan -1 


(1 -a) 2 n a k x s j-k 


(1 - a) 2 cl^x c i — k 
k = 0 


Taking the sine and cosine of this estimate, namely, 65 


y c — cos A a>T s 


where A6kA6k-\ is the kth data symbol phase before 
differential encoding at the transmitter, Aco is the unde- 
sired Doppler term, and N r c k and N ** are input noise 
terms. The noise free sample value signals at the output 
of switches 53 and 54 are respectively applied as one 
input signal each to multiplier junctions 111 and 116. 
These junctions 111, 116 also receive additional input 
signals as processed by Doppler estimator circuit 100. 

Proper detection of A0k requires remo val of the Dop- 
pler from the signals in equation (11). The term Ts/2 is 
a tangent ratio, as is seen by an examination of equation 
(6), and thus we use an arc tangent circuit 95 of any 
well-known type connected to the outputs of digital 
filters 93, 94. Multiplier circuit 96 multiplies the input 
thereto by a factor of 2 to produce the estimate AcbT s 
rather than Ao)T 5 /2. The quadrature Doppler estimate 
signals of equation (8) are cross-correlated by circuits 
110, 117 and 115, 112. These cross-correlated signals are 
processed in the form: 
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2ck = x c k'y c + x s k‘y s 

— P cos [(Aw — A oi)T s + A0*] + N c k cos Aw 7^ + 

N s k sin Awl; 

Zsk = Xsk'y c - Xck'ys 

= P sin [(Aw — Aw) 7^ + A0*] + N s k’ cos Aw T s — 

A/**' sin Aw 7^ 

Finally, the decision variable is determined by 
weighting z c * and z^ with the quadrature coordinate 
projections of any of the possible data symbol phases 
and summing the result to produce 
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common to all baseband and IF equivalent circuitry. 
FIG. 4 requires no further discussion because it uses the 
features and principles described in connection with 

FIG. 3 hereinbefore. 

5 

PERFORMANCE OF THE DOPPLER 
ESTIMATOR 

To assess the performance of the Doppler estimator 
IQ 100 operating in accordance with this invention as de- 
fined in equation (7), we must characterize the statistical 
properties of the noise samples N c £and N s *. Using terms 
(1) and (2) in equation (4), it is straightforward to show 
that 


15 


V; = Z c k COS * 


2777 

M 


• 2777 

+ Zsk sin 




= P cos I (Aw — Aw) 7"* -j- 




f air. + -Vrf' sin f 


Nck cos 


i = 0,1, . . . M - 1 


20 


25 


which then goes to a decision device. If the modulation 
is coded, then the Doppler-corrected signals l c and Q c 
would be input to a maximum-likelihood decoder over 
leads 123, 124. Multipliers 125 and 130 receive the co- 
sine and sine of the various possible phase angles that 30 
are available in the coding format. These reference 
phasors for the decision making process are applied to 
leads 121 and 122 connected to an input of multipliers 
125 and 130. Component 135 is a summer for summing 
the two outputs of multipliers 125 and 130. Summer 135 ^ 
is connected to a decision device as is shown which 
selects the phase angle that is closest to the various 
reference angles. 

FIG. 4 is a baseband implementation of the invention. 
The primary differences between the IF and baseband 
implementations (FIGS. 3 and 4, respectively with like 
numerals designating like components and functions) is 
described as follows. In the IF operation signals are 
delayed and multiplied (31, 32, 28 and 41, 42, 29) at the 
carrier frequency producing sum and difference terms. 
The difference terms represent the desired signals 45 
which contain either the Doppler alone (Eqs. (5) and 
(6)) or the sum of the data and Doppler (Eq. (9)). The 
sum terms represent second harmonics of the carrier 
and are removed by the zonal filters 34 and 36 in FIG. 

3. In the baseband operation, the carrier is first removed 50 
from the input by the in-phase (I) and quadrature (Q) 
demodulators 410 and 411 in FIG. 4 which demodula- 
tors receive the sine and cosine terms of a carrier in 
conventional operation. The demodulation once again 
produces sum and difference terms. The sum terms 55 
(appearing at leads 114 and 115 at the output of filters 
412, 413) are processed using delay and multiply opera- 
tions to produce the desired signals needed for data 
detection and Doppler estimation. One difference, how- 
ever, is that here, in the baseband configuration, we 60 
must sum (at 416, 417, 418 and 419) the I x I rand I x Qr 
signals with the Q x Qrand Q x It, signals respectively, 
to yield the proper signals. (The subscript “T” denotes 
a Trsec or T 5 /2-sec. delayed version of the respective 
signal as appropriate). In the IF configuration, the re- 65 
quired signals are formed from the I x Irand I x Qr 
signals alone. This difference is a consequence of multi- 
plying baseband signals rather than IF signals and is 


N c k = A c k cos Aw7V2 — A s k sin Aw TJ2 
N s k = A c k sin Aw7y2 + A s k cos AwT^/2 


(14) 


where A c k and A s k are, respectively, samples (at 
U=[(4k-f 3)/4]T 5 ) of A c (t) and A^t) defined by 

\l — 

Ac(0 — N P ndt) sin 6(t — T s / 2) + 

M7 rts(t) cos 6(t - Ts/ 2) + n</t — Ts/2) sin 0(0 + 

M"? ris(t — 7V2) cos 0(0 + nAt)ndt - T s /2 ) + n^n^t — T s /2) 
and 

\l — \l — (15b) 

As(t) = N P rtf/t) cos 0(f — Ts/2) — N P rts(t) sin 9{t — T s /2 ) — 

n</t - Ts/2) cos 0(0 + M7 n£t — Ts/2) sin 0(0 + 

n</l)nlt - Ts/2) - n&M - Ts/2) 

The non-Gaussian orthogonal components A c k and A s k 
are uncorrelated. Furthermore, B^ and A c k both have 
zero mean and variance 


cr 2 = PNqB 


( 1 Ab 

l 1 + — J 


06) 


where B is the input IF filter bandwidth. 

Using (14) in (6) and applying straightforward trigo- 
nometric operations, we can write (7) in the normalized 
form 


( 17 ) 


Aw T s — Aa)Ts = 2 tan" 


NIT 


P + z 


= A® 


where 


z = N 1 - a 2 f n a k - ■ £{£} = 0, of = 1 


*=0 <r 


i = N 1 - a 2 2 a k -A±-L ■ £{{} = o, of = 1 
k = 0 a 
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-continued 


PKa ( 

' PK a \ 

!”i + ^£.1 

1 + a 

2<r 2 V 

. ) 

L 1 + ip J 

' Ka “ 1 - a 


5 


10 

(no Doppler) MDPSK, the average symbol error prob- 
ability is given by 

P s = (21) 


Thus, Ao)T s is an unbiased estimator of Aa)Tj. 

Making the practical assumption that K a is large 
actually a value of K a =10 is sufficient) and applying 
central limit therorem-type arguments, then £and £are 
approximately Gaussian distributed. Under this condi- 1° 
tion, the probability density function of 
vAA<I>/2=(Aa) — Aa>)Tr/2 is given by 


Pv(v) = 


e*P (-p) 
27 r 


+ 


(19) 


15 



where E 5 =(log 2 M)E^is the energy per symbol with E b 
the bit energy. 

When Doppler is present, which results in a noisy 
phase A<t>=(Ao> — Ao>)Tj being added to the decision 
variable [see (13)], then using Eqs. (11) and (43) of the 
above-noted Pawula et al. article, equation (21) is modi- 
fied to 


Pst &< ») = 



exp (— p sin 2 v) cos v [1 + erfi \fp" cos v)]; |v| ~ tt 

The variance of the estimator can be computed from the 
second moment of (19) which for K a large becomes 35 


which for AO = 0 reduces to (21). Averaging (22) over 
the probability density function of AO gives the average 
symbol error probability performance of the Doppler- 
corrected differential detector, namely, 



^Acj Ts 



( 20 ) 


... f" 


40 


P£A<t>)p&<p(A<t>)dA<P 


where P y(y) & given by (19) 


-jf 


(23) 

P£2v)p v (v)dv 


This same result can be obtained directly by following 
an approach analogous to that taken in A. J. Viterbi, 
and A. M. Viterbi, “Nonlinear Estimation of PSK- 45 
Modulated Carrier Phase with Application to Burst 
Digital Transmission,” IEEE Transactions on Informa- 
tion Theory, Vol. IT-29, No. 4, July 1983, pp. 543-551. 

EVALUATION OF THE AVERAGE SYMBOL 

ERROR PROBABILITY 50 

In R. F. Pawula, S. O. Rice, and J. H. Roberts, “Dis- 


EVALUATION OF THE AVERAGE BIT ERROR 
PROBABILITY 

To evaluate the average bit error probability, we 
follow the procedure derived by Lee in “Computation 
of the Bit Error Rate of Coherent Mary PSK with Gray 
Code Bit Mapping,” IEEE Transactions on Communi- 
cations, Vol. COM-34, No. 5, May 1986, pp. 488-490 
for ideal MPSK and MDPSK and extend it to the case 
where a noisy phase is added to the decision variable. In 
particular, we define 


A t < A<t>) - * 


in (a<I> - (2 <+ 1 Jy eXp [ -_ ^"( 1 " cos ( A<t> - ( 2i ' + 1 ) - M - ) cosa )] 


(24) 


477 


da - 


sin ^A<t> — (2/ - 


M 


1 — cos ^A<J> — (2/ + 1) )cos a 

) ff exp [-~fr( 1 - C0S ( A<I> - (2 ‘'- 1)- ^) cosa )] 


4t r 


/ 2L e 

V 


— cos ^A<t> - (2/ — 1) JL ^jcos a 


da 


tribution of the Phase Angle Between Two Vectors 
Perturbed by Gaussian Noise,” IEEE Transactions on 
Communications, Vol. COM-30, No. 8, August 1982, 
pp. 1828-1841, it was shown that for uncoded “ideal” 


which represents the probability that the detector out- 
put phase is in the region [(2i-l) /M, (2i + l) /M] under 
the hypothesis that the input phase A0* = O. Then, the 
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conditional symbol error probability of (22) can be 
expressed as 


M- I 

Ps(A<t>) = 1 

/= I 


= / Flu. 2 v)p v (v)dv 


W,/3) = 


sin (j3 — il/) 


*ir T Es 

f exp L-— 


(1 — cos IP — i/i) cos a) 


A'o v ^ ' 

1 — cos 03 — it) cos a 


. , For any fixed fi, F(i//, /3) is easily evaluated using a 

and assuming a Gray code bit assignment with Ham- G auss-Chebyshev formula. In particular, 
ming weights w m> the conditional bit error probability 
is given by 


^ A<t> > = T^aT m h 


20 A 

(26) j = 


rf rf 

§ J(cosa)da = 2 f 

J-f J° 


J{cosa)da 


Here represents the number of bit errors that occur _ j C0S JL L 

if the detector output phase lies in the region corre- 25 Jo J >> 2 / 

sponding to A m (A<t>). 

Recognizing from (24) that Aa/.,(AO)=A,(- A<I>), Letting z = cost}, then 
and since PA<£> (A<f>) is a symmetric probability density 
function then Aa/-;=A; where the overbar denotes 
averaging over PA<t> (A<t>). Hence, when the same 30 

average is performed on (25) and (26), we get £ 1 ■ 


1 -+- COS 7} 
2 


M/2 

P s = 2 2 A m — Am/ 2 

m— 1 


J - 1 Nl-z2 

- N ( 1+1 


1 r ^/2 

g2 M 


(w m + - w.v//2 ^.v//2 


where N is chosen large enough to give the amount of 
accuracy desired and 


2k = cos ( ^ 2k 2 y 1 ^ 77 \ k ^ I, 2, .... N 

Lee, in the above-cited article, showed that for M=4 \ / 

and M — 8, simple closed form expressions could be 

obtained for the bit error probability of ideal MDPSK. i.e., the zeros of the Nth degree Chebyshev polynomi- 
The analogous expressions for the case of corrected 45 nal T.y(x) — cos(Ncos~ 1 x) < 


Doppler present are 




(M = 4) 


Similarly F(U#) of (31) can also be evaluated using a 
Gauss-Chebyshev formula. First breaking the integral 
(29) on ( — 71, 77 ) into two integrals, one on (— 7T,0) and one 
on (0, 77 ), and then using double angle trigonometric 
50 identities, we arrive at the desired result, namely, 




(A/= 8) 


W) “ 77 l ^ { i( z t) 


where Z / is given by (35) with the index k replaced by 
1, and 


g(z) = Pz (z)[F + w,z) + F-(M] (37) 

with P/z) obtained by replacing cos v with z in (19) and 


±(2z N 1 - z 2 ) cos — (2 z 2 — 1) sin 


F ± (ifj,z) = 


X 
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-continued 


* 

* k= 1 


exp^ — — — ^1 — [(2 z 2 — 1) cos »!/ ± (2z N 1 — z 2 ) sin Ui] \ ^ ^ 


1 - [(2 z 2 - 1) cos vi/ ± (22 N 1 - z 2 ) sin vi/] \ 1 ^ ** 


Illustrated in FIGS. 5 and 6 are the average bit error 
probability performances of 4DPSK and 8DPSK as a 
function of bit energy-to-noise ratio, E$/No, in dB with 10 
K a as a parameter. These numerical results were ob- 
tained from equations (29) and (30), respectively. For 
the purpose of comparison, the bit error probability 
performance for perfect Doppler compensation (i.e., 
K a = oo ) is also illustrated in these figures. It should be 15 
observed that to achieve about 0.5 dB requires a value 
of K on the order of eighty (80). This seemingly high 
value of K comes from the fact that the Doppler esti- 
mate must be multiplied by 2 (see the multiplier 96 in 
FIGS. 3 and 4 ) before using that estimate to remove the 20 
Doppler from the data detection signals. This, unfortu- 
nately, increases the variance of the estimator by a fac- 
tor of 4 (see Eq. (20)). The requirement for multiplying 
the Doppler estimate by a factor of 2 is tied to the fact 
that within a given symbol interval, the change in phase 25 
over half a symbol interval is independent of the data 
symbol itself and contains only the Doppler induced 
phase shift, thus allowing its estimation. 

PRACTICAL CONSIDERATIONS 30 

In bandlimited applications, the rectangular baseband 
pulses assumed thus far should be shaped to limit adja- 
cent channel interference (ACI) and control intersym- 
bol interference (ISI). Since 0(t) is no longer constant 
over a symbol time, then in order to still achieve x c k and 35 
x s k of the form in (6), i.e., I and Q Doppler estimation 
samples which are independent of the data modulation, 
we must impose a requirement on the pulse shaping 
such that 

40 

m + DTs) - m + \)T S ) = 0 (39) 

where 0(t) is the filtered version of 0(t). Equivalently, 
we require a pulse shaping such that 

45 

sin m + l)T s ) = sin 9{(k + \)T S ) (40) 

cos §((k + $)T S ) = cos S((k + i)T s ) 

Assuming that one chooses a pulse shaping to satisfy 50 
(40), then, in general, the samples taken at 
t = (k+ l/2)Tj(i.e., the midpoints of the symbol interval 
for data detection purposes will be subject to ISI. Thus, 
the question may arise: Can we choose a pulse shaping 
that satisfies (39) yet still allows data detection without 55 
any penalty due to ISI? This question has a fortuitous 
answer which depends on applying ingenuity which 
relates to additional features of the invention. 

We shall show by reference to equations depicted and 
described relative to FIG. 11 (described in additional 60 
detail in Appendix A hereinafter) that we can obtain the 
characteristics of a 100% raised cosine system by split- 
ting this overall transmission characteristic into a 100% 
root raised cosine filter 1100 at the transmitter and an 
ideal low pass filter 1110 at the receiver combined with 65 
two summed samples (1120, 1125, 1130, 1135) prior to 
demodulation. THe samples are taken one-half symbol 


times. The ideal low pass filter 1110 followed by two 
samples at T 5 /4 (delayed by 1135) and 3T*/4 supplied to 
a summer 1130 yields an equivalent operation to a 100% 
root raised cosine filter. 

Suppose that we consider an overall (transmitter plus 
receiver) pulse shaping of the 100% excess bandwidth 
root raised cosine type. Such pulse shaping is known 
and has a transfer characteristic of the form 


f\a)) = Pj{o>)Pr{ix >) = 


T s cos w7V4; |a)| = 2 n/T s 
0; otherwise 


(41) 


with a corresponding impulse response of 



Since the zero crossings of p(t) occur at 
t= ±(2k+ l)Ty/4; k= 1, 2, 3, . . . , then if p(t) were used 
as the underlying pulse shape in a digital pulse train of 
the form m(t) = 2a w/ ,(t — nT) with {a„} representing a 
multilevel data sequence, then m(t) taken at 
t = (n-f 1/4)T 5 and (n-f-3/4)T* would be ISI free and 
equal in amplitude. Thus, using such an m(t) for the I 
and Q data signals sin 0(t) and cos 6( t) allows (40) to be 
satisfied. 

The question remains as to where to sample and how 
to process the I and Q data signals for data signals for 
data detection purposes. Since, as pointed above, sam- 
ples of sin 0(t) and cos 0(t) taken at t~(n + 1/4)T S and 
(n + 3/4)Tjare ISI free, then it seems reasonable to sum 
these two samples prior to differential detection. Since 
p(t) has value unity at these sample points, i.e., the same 
value as a unit rectangular pulse, then compared to a 
single sample of this rectangular (ISI free) pulse, the 
above sum yields a 6 dB increase in power of the pre- 
detection sample in the absence of noise. 

Ordinarily from the standpoint of maximizing the 
SNR of a single sample of I or Q, one would split the 
total pulse shaping P(ct)) equally between the transmit 
and receive filters with the filters at each location taking 
the form of a square root raised cosine characteristic. 
Note that from an overall system standpoint the filters 
multiply and thus yield an equivalent raised cosine sys- 
tem response. Such a split would result in the noise 
samples at t=(n+ 1/4)T 5 and (n + 3/4)Ts being corre- 
lated. Suppose instead that half of the filtering (square 
root raised cosine) is placed at the transmitter and a 
brick wall filter (ideal low pass) is used at the receiver. 
This results in a 3 dB penalty in the noise power of a 
single sample. However, it can be shown that noise 
samples separated by Ts/2 are now uncorrelated and 
hence the increase in noise power produced by sum- 
ming the above two samples to produce a single pre- 
detection sample is 3 dB. Finally then, for the root 
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the pre-detection sample formed from the sum of sam- 
ples taken at t = (n + 1/4)7^ and (n+3/4)T 5 is identical to 
that of a single sample of a rectangular pulse and thus 
there is no ISI penalty. 

In fact, the above-noted result, which appears rather 
surprising at first, can be supported quite easily by 
showing that, in the absence of Doppler, the split filter- 
ing with an ideal low pass filter at the receiver is indeed 
the optimum filter split for Nyquist-type signaling with 
two nonzero uniform sample points (as is proved by 
Appendix A). In this context, the above scheme bears 
resemblance to the duobinary approach discussed ear- 
lier but is entirely binary and does not depend upon the 
three level complexity of the duobinary approach. 

When Doppler is present, G(a>) of FIG. 12 is replaced 
by 

G(co - Aw) = 1 + exp {-J{o - Aw)7y2} (43) 

s= 1 + a exp {—juTs/2}; a = exp {/Ag>7V/2} 

thus requiring knowledge of the Doppler frequency 
shift. What is used is the Doppler estimate Aa> for Aco in 
the implementation of equation (42) analogous to the 
use of Ao) instead of Ao> for the final Doppler removal. 

If the matched filter follows the differential detector, 
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filter becomes a true matched filter, i.e., a root raised 
cosine filter. Furthermore, the bandwidth of the brick- 
wall filter is narrowed to the symbol rate. Differential 
detector 720 of FIG. 7 operates in the manner described 
5 previously and requires no further discussion. 

A practical implementation of the receiver of FIG. 7 
requires an understanding and use of acquisition 
(“ACQ”) and tracking (“TRK”) modes which will be 
discussed in connection with additional figures such as 
10 FIGS. 7A through 1C. If Aa> = Aco = 0 and the band- 
width of the brickwall filter 740, FIG. 7, is equal to 
1/Tjthen filter 740 is exactly a matched filter for a root 
raised cosine filter with 100% excess bandwidth at the 
transmitter. If Ao> is not equal to zero then filter 740 
15 represents a pseudo-matched filter (e.g., an approxi- 
mately matched filter). The bandwidth of brickwall 
filter 740, in the presence of Doppler frequency shift, 
should be equal to (l/T 5 +f^,MAX). 

FIG. 7A presents additional details of the circuitry 
20 and operation of the Doppler estimator 790 shown sim- 
ply as a block diagram in FIG. 7. A tracking mode for 
the Doppler estimator circuit 790 of FIG. 7A will be 
described. An output from the brickwall filter 740 is 
applied as an input signal to another differential detec- 
3 tor 781 which operates in the manner previously de- 


the circuit operation is once again G(w) of FIG. 12. It 
should be noted that one could indeed also use a differ- 
ential detector followed by the G(<d) combination for 
the data detection portion of the receiver in an effort to 
produce some commonality in structure. However, this 
later approach is suboptimum and leads to a degrada- 
tion in output SNR. Also, in this case, the probability 
density function of the differential detector output 
(which is not Gaussian) is changed by the matched filter 
and hence, the bit error probability performance expres- 
sions of equations (24) through (28) are no longer valid. 

One additional point, alluded to above, needs to be 
discussed. Since in the ACQ mode, the differential de- 
tector operates over a Ty sec. interval, the output of the 
digital filter will proportional to exp {— j Ao>T 5 } rather 
than exp {— j AwTy/2} as in the tracking (TRK) mode 
where the differential detector operates over Ty/2 sec. 
Since implementation of the matched filter for data 
detection requires the complex quantity “a” of equation 
(42), then in the ACQ mode, one must take the square 
root of the digital filter output in order to make this 
quantity available. On the other hand, in the TRK 
mode, where “a” is directly available, then to achieve 
the quantity exp {— j Aa)T<J which is needed for Dop- 
pler removal, we square the digital filter output. 

FIG. 7, described below, relates to differential detec- 
tion with Doppler estimation. This FIG. 7 includes, in a 
practical implementation, the theory discussed above 
concerning data detection and Doppler correction. A 
review of equation 43 shows that if the Doppler amount 
is zero then the matched filter of FIG. 7 is in accor- 
dance with FIG. 12. 


scribed, i.e., that differential detector 781 of FIG. 7A 
differentially detects two samples within each full sym- 
bol interval. The samples are separated by T s /2 in order 
to “wipe off’ the data phase, and thus result in a com- 
plex signal of Doppler induced phase T s /2. The sam- 
pled output from detector 781 is applied to a digital low 
pass filter 782. THis digital low pass filter is comparable 
to those described earlier in connection with FIG. 3 and 
3 5 is employed to improve the estimate of the Doppler by 
removing the effect of noise on Ao)Ts/2. An amplitude 
normalizer 783 normalizes the envelope of digital filter 
output 782 to unity. This normalized output from 783 is 
provided both to the pseudo-matched filter 710, FIG. 7, 
40 and to a squaring circuit 784. A squared output from 784 
is applied to the residual Doppler corrector 730, FIG. 7. 
As previously explained, the residual Doppler is re- 
moved at the output of 730 and the output signal on lead 
731 contains only the differential detected phase data 
45 which is recovered in standard operation by differential 
detector 720, FIG. 7. 

In many applications, Doppler estimation, such as 
that achieved by estimator 790, first takes place via a 
coarse Doppler estimate which is obtained in the pres- 
50 ence of a known data sequence, e.g., all zeros. FIG. 7B 
depicts a fast acquisition circuit that is used when the 
data sequence is known to the receiver. Such preor- 
dained knowledge at the receiver is accomplished by 
having the transmitter send a constant phase for a short 
55 time in the order of about twenty symbols as a preamble 
to every data packet. Such a preamble results in a fast 
acquisition of the Doppler estimate. After that fast ac- 


In FIG. 7, the box legend “C” represents a conjugate 
operation. Pseudo-matched filter 710 includes brickwall 
filter 740, multiplier 750, a half symbol delay 760 and a 
summer 770. Multiplier 750 requires two signals, (1) a 
filtered output signal from brickwall filter 740 and (2) an 
input signal in the form of the complex quantity 
a=exp(j Ao)T 5 /2). The complex quantity exp(j 
AcoTy/2), obtained by Doppler estimator 790, is re- 
quired for Doppler removal in the pseudo-matched 
filter prior to differential detection in detector 720. In 
the absence of Doppler, a= 1 and the pseudo-matched 


quisition the receiver can switch to the tracking mode 
for open loop Doppler estimation where the data phase 
is not known by the receiver. 

In the above-noted acquisition (ACQ) mode, the 
change in 0(t) is known (for the all zeros sequence, this 
change is indeed zero) and thus one would expect that 
the differential detector output samples could be used as 
inputs to a digital filter in the Doppler estimator. Unfor- 
tunately, as we have just pointed out above, the differ- 
ential detector output samples are themselves functions 
of the Doppler estimate and thus implementing the 
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system in the fashion would result in a feedback or 
closed loop system. Since one of the goals of our inven- 
tion was to maintain a feedforward (open loop) configu- 
ration, we compromise a bit by reversing the order of 
matched filter and differential detector operations for 5 
the Doppler estimation portion of the receiver in the 
ACQ mode. 

Again in FIG. 7B, the brickwall filter 740 is used to 
noise limit and has a bandwidth that is equal to the 
symbol rate is Hertz plus the maximum frequency shift 10 
that is to be expected. A pair of differential detectors 
721 and 722, operating as previously described, include 
switches 727 and 728 to obtain two samples during each 
symbol interval. THe noise in these samples at 1/4T* 
and 3/4Tyis independent of each other but both include 15 
the same Doppler-induced phase. Accordingly, adding 
these samples in summer 729 enhances the signal-to- 
noise ratio by approximately 3 dB as explained earlier. 
The complex signal at the output of summer 729 in- 
cludes a noisy version of the Doppler-induced phase 20 
and the version is passed through a digital filter 782 and 
a normalizer 783. Two output signals are developed. 
One output signal is the exp(jA<uTj) term and the other 
output signal is fed through a square root circuit 734, 
producing exp(jAo>T 5 /2). 25 

FIG. 1C is an equivalent fast acquisition circuit 
which employs only one detector 742, a delay and add 
circuit 743, (including sampler 729), digital filter 782, 
normalizer 783 and square root circuit 734. The opera- 
tion of the circuitry of FIG. 7C is functionally equiva- 30 
lent, in fast acquisition, to that of FIG. 7B and thus 
requires no further explanation. 

Only a brief correlation of the circuit components of 
FIG. 8 relative to the circuits of FIGS. 7, 7A, 7B or 7C 
need be described. FIG. 8, a simple block diagram of a 35 
dual mode receiver in complex form, is essentially the 
integration of the circuits shown in FIGS. 7A and 7C. 

In FIG. 8 some circuit elements are shared between the 
circuits of FIGS. 7A and 7C to reduce the complexity 
of the dual mode receiver of FIG. 8. Switches 860, 873 40 
and 875 are switches to select the acquisition mode or 
the tracking mode. Only one mode at a time is em- 
ployed. FIG. 8 components 805, 810, 820, 815 represent 
the pseudo-matched filter 710 of FIG. 7. Components 
880, 881, 890 and 891 represent the data differential 45 
detector 720 of FIG. 7. Components 892, 893 represent 
the residual Doppler correction circuit 730 of FIG. 7. 
Components 810, 826, 850, 855, 782, 783, 784 represent 
the Doppler estimator in tracking mode of FIG. 7A. 
Components 805, 810, 840, 845, 830, 835, 836, 838, 782, 50 
783, 734 represent the fast acquisition circuit of FIG. 

1C. The operation of FIG. 8 requires no further descrip- 
tion. 

FIG. 9 is an actual baseband I-Q implementation of 
the Doppler-corrected differential detector analogous 55 
to FIG. 2 but including all of the above practical con- 
siderations for the bandlimited application as incorpo- 
rated in the simple illustration of FIG. 8. As described 
earlier, the elements 31, 32, and 51 are each delays of 
one-half symbol interval as are elements 41, 42 and 52. 60 
The eye diagram (FIG. 10) of the brickwall filter out- 
puts corresponding to sin 0(t) and cos B(t) is illustrated 
for 8DPSK. Some other practical simplifications have 
been made in regard to the manner in which the tan -1 
operation is performed. In particular, the quadrature 65 
signals necessary to remove the Doppler (see FIG. 9) 
can be obtained by noting the following trigonometric 
identities. 


'f )= UB 


Also, allowance has been made for the possibility of 
acquisition and tracking modes of operation as dis- 
cussed above. 

SIMULATIONS RESULTS 

A software simulation of the Doppler frequency esti- 
mator was conducted and proved: 1) verification of the 
analytical results, and 2) inclusion of the effects of inter- 
symbol interference (ISI) due to timing jitter on the 
performance of the system under consideration. The 
simulation and analytical results were in excellent 
agreement as shown by FIGS. 13 and 14 included 
herein. 

APPENDIX A 

On the Optimal Filter Split for Nyquist-Type Signaling 
with Two Nonzero Sample points 

In the foregoing, it was shown that 100% excess 
bandwidth root raised cosine signaling could be used to 
satisfy the Doppler estimation requirement of equal, 
zero ISI sampling points spaced T s /2 apart and symmet- 
rically placed with respect to the midpoint of the sym- 
bol interval. From the standpoint of optimum data de- 
tection, the conventional Nyquist approach is to use an 
equivalent flat channel (i.e., the sum of all translates of 
the overall frequency response should be a constant 
over the Nyquist bandwidth) and comparison of a single 
sample to a threshold. Raised cosine signaling achieves 
such an equivalent Nyquist channel with a single ISI- 
free sampling point. Thus, to use this technique in com- 
bination with the above Doppler estimation circuit, one 
would require an additional root raised cosine filter at 
the receiver. (Actually, for the baseband implementa- 
tion, two such additional filters would be required; one 
for the I channel and one for the Q channel). 

To avoid the necessity of this additional root raised 
cosine filter, the overall channel response remains root 
raised cosine, as required for the Doppler estimate, and 
for data detection the pair of samples taken per symbol 
are algebraically summed and used as a single sample 
for comparison with a threshold. Using the above filter 
split (i.e., all of the filtering at the transmitter and none 
at the receiver), it was shown that the error probability 
performance achieved was indeed equal to that using 
the conventional Nyquist approach (i.e., raised cosine 
signaling with single sample per symbol detection) as 
described above. 

In this appendix, by interpreting the summation of 
two Ts/2-spaced samples as a simple digital filter, we 
show that this choice of filter split is indeed optimum in 
accordance with the traditional notion of putting equal 
amounts of the overall channel response in the transmit- 
ter and receiver. 

SYSTEM MODEL 

Consider the mathematical model of the system illus- 
trated in FIG. 11 . Actually, this system represents only 
the in-phase (I) or quadrature channel (Q) of the base- 
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band receiver implementation. For our purpose here, it 
is sufficient to consider only either the I or Q channel. 
The data source represents either the I or Q component 
of the MPSK modulation and as' such produces M/2- 
level symbols at a rate l/T s . These symbols are trans- 5 
mitted (we consider only the baseband portion of the 
overall system) as a pulse sequence mT(t) with pulse 
shape p7<t) equal to the inverse Fourier transform of the 
transmit filter Pz<aj). At the receiver the noise cor- 
rupted sequence x(t)=m7<t)-hn(t) is passed through the 
receiver filter P/j(o)) to produce y(t) = m(t) + n(t). As 
discussed above, it is desired that m(t) be a pulse se- 
quence with a pulse shape p(t) such that samples of m(t) 
taken at t = (n+ l/4)T 5 and t = (n+3/4)Ty result in zero 
ISI. To achieve this, we have selected P(oi) and p(t) as 15 
in equations (41) and (42), respectively. 

The receive filter output y(t) is sampled twice per 
symbol as shown in FIG. 11 and the sum of these two 
samples is inputted to the differential detector. From 
the standpoint of analysis, the two samplers (A, 1110 20 
and B, 1125) and the T 5 /2 delay 1135 can be replaced by 
a simple digital filter G(a>) 1200 and a single sampler 
1215 as shown in FIG. 12. Thus, by analogy with a true 
Nyquist channel, wherein a single sample of the output 
is used for purposes of data detection, we can regard 25 
samples of z(t) (see FIG. 12) as the pre-detection vari- 
ables and the cascade of P j r(cj) and the digital filter 
G(ca) as the total receive filter. 

It is well known that for additive white Gaussian 
noise (AWGN) channels with no ISI, i.e., Nyquist pulse 30 
shaping, the average symbol error probability is mini- 
mized (with respect to the noise) by selecting identical 
transmit and receive filters each of whose transfer func- 
tion has a magnitude response equal to the square root 
of the overall transmission characteristic. Thus, in terms 35 
of the equivalent representation of FIG. 11, optimum 
performance would be achieved by selecting 
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claims when read in light of the description and draw- 
ings. 

What is claimed is: 

1. A receiver in a communication system which pro- 
duces a Doppler effect in an incoming multiple differen- 
tial phase-shift-keyed signal having a given symbol in- 
terval, with the receiver for the incoming signal having 
a differential detector for restoring the data by multiply- 
ing the incoming signal with a full symbol interval de- 

10 layed version of itself (or an equivalent sampled version 
thereof) and applying the restored signal to the detec- 
tor’s output, said receiver comprising: 

means for delaying the incoming data signal (which 
has the Doppler effect therein) by a half symbol 
interval; 

means in an open-loop forward acting circuit for 
estimating the amount of the Doppler effect from 
the incoming data signal and the half-symbol de- 
layed signal; and 

means for subtracting the estimated Doppler amount 
from the detector’s output, 

2. A communication system having a transmitter and 
the receiver of claim 1 separated by a bandlimited me- 
dium over which the transmitter sends to the receiver, 
a data modulated multiple differential phase-shift-keyed 
signal having a predetermined symbol interval, T s , con- 
taminated by said Doppler effect, said system compris- 
ing: 

signal shaping split between the transmitter and the 
receiver to yield an impulse response at the re- 
ceiver in the form of 



|/> 7 <a>)| 2 = j/Mo>)| 2 |G(m)| 2 = |Ao>)| 2 /|/>/Ka>)j 2 


(A-l) 


However, the magnitude squared transfer function of 
the digital filter G(o>) is given by 
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|G(a>)| 2 = (2 cos W7V/4) 2 for all oj 


(A-2) 


45 


Then using (41) and (A-2) in (A-l) gives the desired 
result, namely, 


I Pr(&) I 


- 0 ^. 


(A-3) 


M ^ lir/Ts 


0; otherwise 
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|w-)| f cos o)7y4; [to | = 2tt/T s 

l 0; otherwise 


The above description presents the best mode con- 
templated in carrying out our invention. Our invention 60 
is, however, susceptible to modifications and alternate 
constructions from the embodiments shown in the 
drawings and described above. Consequently, it is not 
the intention to limit the invention to the particular 
embodiments disclosed. On the contrary, the invention 65 
is intended and shall cover all modifications, sizes and 
alternate constructions falling within the spirit and 
scope of the invention, as expressed in the appended 


where t is time and T s is the symbol duration. 

3. A communication system in accordance with claim 

2 wherein said differential detector for decoding the 
data represented in the received signal has the Doppler 
effect in the detector’s output, said system further com- 
prising: 

means in the detector including said receiver’s signal 
shaping for sampling the received signal at times 
t — (n + 1/4) Tyand (n + 3/4) T^to obtain intersym- 
bol interference-free (“ISI”) data signal samples, 
wherein “n” corresponds to the nth symbol inter- 
val of duration T s and n is any arbitrarily selected 
integer. 

4. A communication system in accordance with claim 

3 wherein said means for removing the Doppler effect 
from the differential detector output, further comprises: 

means for multiplying said ISI samples with the Dop- 
pler estimating amount from said subtracting 
means as a Doppler-correcting error to cancel the 
Doppler effect from the differential detector’s out- 
put. 

5. A communication system in accordance with claim 
1 wherein the Doppler effect introduces a Doppler- 
induced phase change that tends to mask differential 
phase changes encoded in said incoming data signal, and 
said system includes a receiver which is further charac- 
terized in that: 

said Doppler effect estimating means is responsive to 
the change in phase of said incoming signal over a 
full symbol interval, which incoming signal con- 
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tains both the differentially encoded data and the 
Doppler-induced phase change, whereas the same 
phase change over half a symbol interval contains 
only the Doppler-induced phase change. 

6. A communication system in accordance with claim 5 
2 wherein said signal shaping is provided by transmit 
and receive filters with the transmitter filter taking the 
form of a square root raised cosine characteristic, the 
improvement in said system comprising: 

a filter at the receiver having a characteristic 
matched to the transmitter’s filter and also exhibit- 
ing a square root raised cosine characteristic that is 
formed from an ideal low pass filter with its output 
connected in tandem to a summer, which summer 
receives the low pass filter’s output, and a half 15 
symbol delayed version of that filter’s output. 

7. A communication system in accordance with claim 
6 wherein the symbol interval is expressed as T*and the 
transmitter’s filter has 100% excess bandwidth and fur- 

2Q 

ther wherein said receiver’s 

ideal low pass filter has a bandwidth equal to l/T s 
expressed in Hertz. 

8. A communication system in accordance with claim 
6 wherein the symbol interval is expressed as T s and ' 
wherein the transmitted signal experiences a maximum 
Doppler frequency shift (“Umax*') and further wherein 
said receiver’s 

ideal low pass filter has a bandwidth of l/T s (in H 2 ) 

plus Umax- 30 

9. A communication system in accordance with claim 
8 wherein the low pass filter’s output is connected to a 
Doppler correcting circuit, and wherein said system is 
further characterized: 

said means in an open loop forward acting circuit 35 
which receives the low pass filter’s output, for 
estimating said Doppler-correcting error signal, 
which error signal is indicative of the amount of 
Doppler correction that is required at the detec- 
tor’s output. 40 

10 . A communication system in accordance with 

claim 9 wherein the data at the transmitter is encoded as 
multiple differential phase-shift-keyed changes in each 
of said symbol intervals, and further wherein said re- 
ceiver includes a data detector; 45 

said data detector including a full symbol delay T s 
and a multiplier circuit for multiplying the incom- 
ing data signal with a full-symbol-interval delayed 
version of itself. 

11 . A communication system in accordance with 50 
claim 9 wherein the Doppler correcting circuit further 
comprises, as said half symbol delay means, a: 

half symbol delay circuit for delaying the incoming 
data signal by a delay time of one-half (TV2) of said 
symbol interval; 55 

means for multiplying that half delayed signal with 
the incoming data signal; and 

means for further processing the multiplied signal 
from said multiplying means before said signal is 
applied to said detector’s output as a Doppler error 60 
correction signal. 

12 . The communication system of claim 11 wherein 
the Doppler error processing means, in said Doppler 
correcting circuit of claim 11 , further comprises: 

means for sampling the received signal at times 65 
t=(n+ + 1/4) T 5 and (n + 3/4) T Jf where “n” corre- 
sponds to the nth symbol interval of duration T s 
and n is any arbitrarily selected integer, and 
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means for multiplying the sampled signals from said 
sampling means with said detector’s output. 

13 . A communication system having a transmitter and 
a receiver separated by a bandlimited channel over 
which the transmitter sends to the receiver, a multiple 
differential phase-shift-keyed signal (“MDPSK”) hav- 
ing a predetermined symbol interval, T*, contaminated 
by a Doppler effect, said system comprising: 

signal shaping means split between the transmitter 
and the receiver to yield an impulse response at the 
receiver in the form of 



where t is time and T s is the symbol duration, 
an ideal low pass filter in said signal shaping means 
for filtering the incoming signal, 
means in said receiver for sampling the received and 
filtered incoming signal at times t = (k-j-l/4) T s 
where k is equal to 1, 2, 3 etc. and (k + 3/4) T s to 
obtain Doppler error samples, and 
an open loop feedforward circuit in said receiver 
connected to said sampling means for removing the 
Doppler effect from the incoming signal in accor- 
dance with said Doppler error samples from said 
sampling means. 

14 . A system in accordance with claim 13 wherein the 
open loop feedforward circuit further comprises: 

means delaying the incoming signal with the Doppler 
effect therein, by a half symbol interval for com- 
parison with the incoming signal itself to form the 
Doppler removal signal. 

15 . An improved communication system of the type 
that experiences a Doppler effect in a signal which is 
being transmitted over a link between the system’s mod- 
ulator and demodulator, said system including a re- 
ceiver with an input for the Doppler-affected incoming 
signal and an output terminal, for supplying decoded 
data, and comprising: 

a modulator at the system’s transmitter for impressing 
information as differential phase-shift-keyed 
changes on a transmitted carrier as an output sig- 
nal; 

a Doppler-inducing communication medium that 
causes a Doppler phase error to become embedded 
in the transmitted signal; 

a system receiver operative on the incoming signal 
for both removing the Doppler phase error and 
restoring the information from the differential pha- 
ses in said incoming signal, said receiver including: 
a demodulator connected to receive the incoming 
signal and restore therefrom the transmitted infor- 
mation which still has the Doppler phase error in 
the restored signal; 

a sensing circuit receiving the incoming signal and 
responsive to a particular portion of each symbol 
interval of the incoming signal having only the 
Doppler-phase therein for deriving a Doppler 
error correction signal; and 
a forward acting open loop correction circuit con- 
nected to said sensing circuit for removing the 
Doppler phase error from said restored signal so 
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that information only is available at the receiver’s 
output terminal. 

16. An improved communication system for commu- 

nicating information, encoded by differential phase 
changes in symbol intervals, over a communication link 5 
which introduces a Doppler-induced phase change that 
tends to mask the encoded differential phase changes, 
said system improvement being based upon the recogni- 
tion that, at a receiver, the change in phase of a received 
incoming signal over a full symbol interval contains 10 
both the differentially encoded data and the Doppler- 
induced phase change, whereas the same phase change - 
over half a symbol interval contains only the Doppler- 
induced phase change, and said system improvement 
comprising at the system’s receiver: 15 

an input signal receiving terminal; 

means connected to said terminal and operative over 
a half symbol interval in response to said input 
signal for deriving therefrom a Doppler estimate 
signal that is indicative of the Doppler-induced 20 
phase change which is present in the incoming 
signal; 

means connected to said input terminal and operative 
over a full symbol interval in response to said input ^ 
signal for deriving therefrom a detected output L 
signal which is indicative of the data encoded 
phase change and including the Doppler-induced 
phase change as an undesired part of the detected 
output signal; and ^ 

means applying the Doppler estimate to the detected 
output signal for canceling the undesired Doppler- 
induced phase change from the detected output 
signal. 

17. A communication system for symbol interval data 35 
in an environment with Doppler frequency shift and 
having a signal transmitter and a signal receiver con- 
nected together via a communication medium, and hav- 
ing transmit and receive filters with the transmitter 
filter taking the form of a square root raised cosine 40 
characteristic, the improvement in said system having 
an estimate of said Doppler frequency shift and com- 
prising: 

a filter at the receiver having a characteristic pseudo 
matched to the transmitter’s filter and also exhibit- 45 
ing a square root raised cosine characteristic only 
in absence of Doppler shift, said filter being formed 
from an ideal low pass filter with its output con- 
nected in tandem to a summer, which summer 
receives the low pass filter’s output, and a half 50 
symbol delayed version of that filter’s output, 
which half symbol delayed version is corrected by 
said Doppler estimate. 

18. An improved communication system of claim 17 

and further comprising: 55 

means at the transmitter for transmitting, as part of a 
data signal, a known sequence of phase shifts; and 

means at the receiver responsive to that known se- 
quence for acquiring a Doppler error signal during 
the time of receipt of that known sequence. 60 

19. An improved communication system of claim 18 
and further comprising: 

means at the receiver responsive to said Doppler 
error acquiring means for condition achieving 
Doppler correction during the receipt of unknown 65 
data phase shifts by said receiver. 

20. A communication system in accordance with 
claim 17 wherein the symbol interval is expressed as T s 
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and the transmitter’s filter has 100% excess bandwidth 
and further wherein said receiver’s: 

ideal low pass filter has a bandwidth equal to 1/T 5 
(expressed in H z ). 

21. A communication system in accordance with 
claim 17 wherein the symbol interval is expressed as T s 
and wherein the transmitted signal experiences a maxi- 
mum Doppler frequency shift (“f d,MAX”) and further 
wherein said receiver’s: 

ideal low pass filter has a bandwidth of lT 5 (in H z per 
second) plus f d,MAX> 

22. A communication system in accordance with 
claim 21 wherein the low pass filter’s output is con- 
nected to a Doppler correcting circuit, and further 
wherein said Doppler correcting circuit is character- 
ized by: 

means in an open loop forward acting circuit which 
receives the low pass filter’s output for estimating 
therefrom, a Doppler error signal, which error 
signal is indicative of the amount of Doppler cor- 
rection that is required; 

a data detector for restoring data at the receiver and 
making the restored data available at the detector’s 
output; and 

means for removing the Doppler from the restored 
data of said detector by applying said Doppler 
error signal from said estimating means to the de- 
tector’s output. 

23. A communication system in accordance with 
claim 22 wherein the data is encoded as multiple differ- 
ential phase changes in said symbol intervals, and fur- 
ther wherein: 

said data detector is a multiple differential phase de- 
coder for restoring the data in a data output signal 
that is formed by multiplying, at a multiplier cir- 
cuit, the incoming data signal with a full symbol 
interval delayed version of itself. 

24. A communication system in accordance with 
claim 23 wherein the Doppler correcting circuit further 
comprises: 

means for delaying the incoming data signal by a 
delay time of one-half of the symbol interval’s 
means for multiplying that half delayed signal with 
the incoming data signal; and 

means for further processing the multiplied signal 
from said multiplying means before said signal is 
applied to said detector’s output as a Doppler error 
correction signal. 

25. The communication system of claim 23 wherein 
the processing means further comprises: 

means for sampling the multiplied signal at times 
t = (n H — 1- l/4)T 5 and (n+3/4) T s , where “n” corre- 
sponds to the nth symbol interval of duration T* 
and n is any arbitrarily selected integer; and 

means for multiplying the sampled signals from said 
sampling means with said detector’s output. 

26. A receiver in a communication system which 
produces a Doppler effect in an incoming multiple dif- 
ferential phase-shift-keyed signal having a given symbol 
interval, with the receiver for the incoming signal hav- 
ing a differential detector for restoring the data by mul- 
tiplying the incoming signal with a full symbol interval 
delayed version of itself (or an equivalent sampled ver- 
sion thereof), said receiver comprising: 

means for delaying the incoming data signal (which 
has an unwanted Doppler effect therein) by a half 
symbol interval; and 



5,007,068 


25 

an open-loop forward acting circuit means connected 
to said delaying means for deriving an estimated 
amount of the Doppler effect from the signal to be 
restored, the estimated amount then being used to 5 
cancel the unwanted Doppler. 

27. A communication system having a transmitter and 
the receiver separated by a bandlimited medium over 
which the transmitter sends to the receiver, a data mod- 
ulated multiple differential phase-shift-keyed signal or a 10 
multiple phase-shifted-keyed signal, with each signal 
having a predetermined symbol interval, T s , contami- 
nated by an unwanted Doppler effect that may be re- 
moved by samples, a half symbol-interval apart, from 15 
the incoming signal, said system comprising: 

signal shaping split between the transmitter and the 
receiver to yield an impulse response at the re- 
ceiver, where said signal shaping is connected, 
with said impulse being in the form of 



where t is time and T*is the symbol duration, and said 30 
signal shaping yields samples for Doppler removal that 
are a half symbol duration apart; and 
means connected to said signal shaping means and 
responsive to said samples yielded thereby for cor- ^ 
recting for said unwanted Doppler effect. 

28. The communication system of claim 27 wherein 
the Doppler error is removed by samples taken a half 
symbol duration apart for Doppler correction, said 
system further comprising: 

means for sampling the multiplied signal at times 
t = (n-f* 1/4)T S and (n + 3/4)T 5 , where n corre- 
sponds to the nth symbol interval of duration T* 
and n is any arbitrarily selected integer, in order to 45 
obtain two independent samples; and 
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Doppler correcting means responsive to said inde- 
pendent samples for correcting for said unwanted 
Doppler. 

29. A receiver in a communication system which 
produces a Doppler effect in an incoming multiple 
phase-shift-keyed signal having a given symbol interval, 
with the receiver for the incoming signal having a 
matched filter detector for restoring the data after Dop- 
pler removal, said receiver comprising: 

means for delaying the incoming data signal (which 
has the Doppler effect therein) by a half symbol 
interval; and 

an open-loop forward acting circuit means connected 
to said delaying means for estimating the amount of 
the Doppler effect from the signal to be restored, 
the estimate of which is then used to cancel the 
Doppler. 

30. A communication system for symbol interval data 
having a signal transmitter and a signal receiver con- 

20 nected together via a communication medium, and hav- 
ing transmit and receive filters with the transmitter 
filter taking the form of a square root raised cosine 
characteristic, the improvement in said system compris- 
ing: 

a filter at the receiver having a characteristic 
matched to the transmitter’s filter and also exhibit- 
ing a square root raised cosine characteristic that is 
formed from an ideal low pass filter with its output 
connected in tandem to a summer, which summer 
receives the low pass filter’s output, and a half 
symbol delayed version of that filter’s output, 
means at the transmitter for transmitting, as part of a 
data signal, a known sequence of phase shifts; and 
means at the receiver responsive to that known se- 
quence for acquiring (“ACQ”) a Doppler error 
signal during the time of receipt of that known 
sequence. 

31. An improved communication system of claim 30 
and further comprising: 

40 means at the receiver responsive to the acqui-sition of 
the Doppler correction by said Doppler error ac- 
quiring means for assuming a tracking condition 
wherein the Doppler correction takes place during 
the receipt of unknown data phase shifts by said 
receiver. 
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